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Abstract. We describe recent theoretical work on the final stage of accretion 
when material passes from an accretion disk onto a white dwarf surface. Our 
calculations focus on understanding the latitudinal spreading and differentially 
rotating profile of this material, so we call it the "spreading layer" (SL) model. 
The SL typically extends to an angle of #sl ~ 0.01 — 0.1 radians, with respect to 
the equator. At low accretion rates (M <J 10 18 g s _1 ) the amount of spreading 
is small, so that the dissipated energy is radiated back into the accretion disk. 
When M is high, such as in dwarf novae, symbiotic binaries, and supersoft 
sources, the material spreads to higher latitudes to be more easily observed. The 
SL may contain shallow surface modes, and we propose that such waves could 
produce dwarf nova oscillations (DNOs). This hypothesis naturally explains 
many key properties of DNOs, including their frequency range, sinusoidal nature, 
sensitivity to M, and why they are only seen during outburst. 



1. Introduction 

If a white dwarf (WD) is weakly magnetized, accreted material reaches its surface 
at the equator and then must lose angular momentum and kinetic energy before 
settling onto the star. When the WD spin is much less than breakup, nearly 
half of the accretion luminosity is released in this process, making this region 
as luminous as the accretion disk and crucial to understanding the luminosity 
from cataclysmic variables (CVs). This transition from accretion disk to WD is 
typically treated in a boundary layer (BL) model. One-dimensional BL studies 
follow the radial flow of material and assume a vertical scale height (e.g. Popham 
& Narayan 1995), while other BL models study both radial and latitudinal 
directions simultaneously using numerical techniques (e.g. Kley 1989a, 1989b). 
Inogamov & Sunyaev (1999; hereafter IS99) approached this problem from a 
new angle to study accreting neutron stars (NSs). Their method follows the 
latitudinal flow of matter and provides information about the spreading area of 
hot, freshly accreted material which is not captured in BL models. We apply 
these same methods to the case of WDs (Piro &: Bildsten 2004a; hereafter PB04a) 
and call this model the spreading layer (SL), to differentiate it from BL studies. 

The SL is found to be a thin hot band of extent #sl ~ 0.01 — 0.1 radians, 
with an effective temperature of ~ (2 — 5) x 10 5 K for M = 10 17 — 10 19 g s" 1 , 
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implying that it contributes to accreting WD spectra in the extreme ultraviolet 
(EUV) and soft X-rays. At low M (<; 10 18 g s -1 ) the spreading is small, so that 
its radiation is absorbed by the accretion disk or its winds. Nevertheless, the SL 
may be important for calculating the underlying continuum scattered by these 
winds. 

The material in the SL is much hotter in temperature and lower in density 
than the underlying WD. This contrast allows waves in the SL to travel freely, 
unencumbered by the material below. We propose that dwarf nova oscillations 
(DNOs) are shallow surface waves in a layer of recently accreted material con- 
fined to the WD equator (Piro & Bildsten 2004b; hereafter PB04b). This model 
explains the main properties of DNOs, including their sinusoidal nature, depen- 
dence on M, large pulsed amplitude in the EUV, and why they are only seen 
during outburst. 



2. The Spreading Layer Calculation 

The SL is treated as a one-zone, plane-parallel layer with 6 (the angle from the 
equator) as the independent coordinate. The surface gravitational acceleration 
is significantly decreased by centrifugal effects due to the nearly Keplerian speed 
of the layer, so we use an "effective" gravitational constant 

GM vj v\ GM v\ 

where vg and v<f, are the latitudinal and azimuthal velocities, respectively. Hy- 
drostatic balance integrates to P = g c sy, where y is the column depth (cgs units 
of g cm -2 ). The flux in the layer is assumed constant since most of the viscous 
dissipation occurs at the layer's base. Radiative diffusion integrates to 

where a is the radiation constant and k = 0.34 cm 2 g -1 is the Thomson scatter- 
ing opacity for a solar composition. 

We insert this one- zone model into a set of differential equations describing 
the steady-state conservation of mass, momentum, and energy, as done by IS99 
for NSs. The total viscous stress, r, is from the turbulence created as high 
entropy fluid quickly rotates against the lower entropy WD material below, so 
parametrizing r in terms of a unitless constant, usl, 

t = a S LPV 2 , (3) 

where v 2 = Vq+v"^. We treat «sl as a free parameter and study its effect on the 
spreading properties for the values «sl = 10~ 4 — 10~ 2 . This range is chosen for 
the physically realistic solutions that result (see PB04a for further discussion). 

Integration requires setting boundary conditions for 6, T, vg, and v^. Shakura 
& Sunyaev (1973; for gas pressure dominating radiation pressure and a Kramer's 
opacity) show the disk subtends an angle at the WD surface 

* dlsk = 1.8 x 10-a^Vr Mr 3/8 < 8 (4) 
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where adisk is the viscosity parameter for the accretion disk, adi s k,2 = ^disk/10 2 > 
Mi 7 = M/(10 17 g s" 1 ), Mi = M/M Q , R 9 = R/(W 9 cm), and we set the factor 
[1 — (r/i?) -1 / 2 ] ~ 1. The material which comes in from the disk must start at an 
angle <; #disk> so we use an initial angle of w 10~ 3 . The initial is set to 0.99-uk, 
where wk = (GM/ R) 1 / 2 is the Keplerian velocity, and the initial T is set using 
the disk's midplane temperature (Shakura & Sunyaev 1973). Finally, the initial 
vq is set by the natural tendency for the differential equations to asymptote to 
the limit where initially the viscous dissipation is radiated away locally (PB04a). 
This is different than what IS99 find for NSs, where more advection takes place, 
due to M closer to the Eddington limit. 

In Figure H we plot the velocities, T, T e g and the pressure scale height, 
h = Pj (pg e g), of the SL for «sl = 1CP 3 and a range of M. At the high accretion 
rates (M £ 10 19 g s _1 ) more advection takes place, resulting in solutions which 
look similar to what IS99 found for NSs. Most notably, the T profiles show large 
peaks at high latitudes due to energy being advected away from the equator 
and deposited higher up on the star. The panel plotting T e g has an additional 
dashed line at 3 x 10 4 K, a fiducial temperature of the underlying accreting WD 
(Townsley & Bildsten 2003) to show the contrast between it and the SL. 

Using the results from Figure 1, we plot the dependence of #sl on M in 
Figure 2. This shows that high M are needed for any hope of appreciable 
spreading to take place, so that typically the SL will be absorbed and re-radiated 
by the accretion disk and its winds. Nevertheless, even when the accretion disk 
covers the main portion of the SL, the profiles of T e g shown in Figure 1 are 
important for understanding the WD spectra. 



3. Comparisons to Dwarf Novae in Outburst 

The high M needed for appreciable spreading suggests that a good candidate to 
show a SL are dwarf novae in outburst. Evidence for spreading can be explored 
by following T e g versus M in the EUV. Since we find that the radiating area 
increases with M, the SL model predicts a shallower change in T e g with M 
than if the radiating area was fixed. PB04a use this relation in an attempt test 
for evidence of spreading during an outburst of SS Cyg (Wheatley, Mauche & 
Mattei 2003), but are not able to make any definite conclusions with the current 
data. Further detailed observations following CVs in outburst are needed to test 
the SL model in better detail. 

The ease of investigating the SL depends on the angle of the observed 
system. Face-on systems show a continuum most likely from the SL region, but 
with strong absorption lines from the disk winds (Mauche 2004). Assuming the 
fractional emitting area is approximately equal to #sl> modeling of this spectra 
finds 6 S l ~ 5.6tf 5 x 10" 3 (for SS Cyg; Mauche 2004), which, though not highly 
constraining, is consistent with our calculations. Edge-on systems show evidence 
of a SL emission that is scattered into the line of sight by line emission in the 
disk winds (OY Car; Mauche & Raymond 2000), but is otherwise obscured by 
the accretion disk. It is therefore difficult to make direct comparisons with the 
SL model. 
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Figure 1. The velocities, temperature, effective temperature and pressure 
scale height of the SL for a WD with M = 0.6M Q and R = 9 x 10 8 cm, 
with a S L = 10~ 3 and a disk = 0.1. From left to right the M are 10 17 , 10 18 , 
10 19 and 3 x 10 19 g s" 1 . The dashed lines shows T off = 3 x 10 4 K, a fiducial 
temperature for the underlying accreting WD (Townsley & Bildsten 2003). 



4. Dwarf Nova Oscillations 



DNOs are oscillations observed during dwarf nova outbursts that have periods 
of P 3 — 40 s, scale monotonically with EUV luminosity, and thus M (Mauche 
1996), and are fairly coherent (Q ~ 10 4 — 10 6 ). Their large pulsed fraction in 
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Figure 2. The spreading angle, 6>sl, versus M for a WD with M = 0.6M Q 
and R = 9 x 10 9 cm and three values of a$L- The spreading angle is defined 
as the angle at which v,p = 0.1z;k- This shows that high M are needed for the 
SL to extend beyond the angle made by a Shakura-Sunyaev disk, equation 
(4) (dashed line). 



the EUV implies that they originate close to the WD surface and have a similar 
covering fraction as the SL (Cordova et al. 1980; Mauche & Robinson 2001; 
Mauche 2004). We propose that DNOs are produced by nonradial oscillations 
in the SL (PB04b). In the following sections we present a simple model to 
describe these oscillations and make comparisons with observations. 

4.1. Nonradial Oscillations in the SL 

Patterson (1981) showed that (at the time) all known DNOs on WDs with 
measured masses have P ;> Pk = 2vr(i? 3 /GM) 1 / 2 . This relation is an important 
constraint for any explanation of DNOs, so we first show why a SL mode should 
mimic such a period. The large entropy contrast between the surface layer 
and the underlying material confines nonradial oscillations to high altitudes 
with little or no pulsational energy extending deeper into the WD. When the 
horizontal wavelength is much greater than the layer depth, these shallow surface 
waves have a frequency 

= 9cf[hk 2 , (5) 

where k is the transverse wavenumber. The SL can be thought of as a waveguide 
with latitudinal width 29q,\ j R and azimuthal length 2ttR, but since #sl <C 1 the 
latitudinal contribution dominates, so k ~ 1/(#sl-R)- Setting g e g = XGM/R 2 , 
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where A <J 1 is a dimensionless parameter that depends on the spin of the layer, 
we rearrange the terms in equation (5) to find 

(GM\ 1 ' 2 ( Xh y /a 

so that the mode's frequency is the Keplerian frequency times a factor less than 
unity (as long as Xh <; 9^ L R), and therefore consistent with Patterson (1981). 

To find how these modes scale with M and M we consider a simple model to 
estimate h. Approximately half of the accretion luminosity is released at the WD 
surface, so that the flux of this layer is given by 47t6> S lP 2 P = GMM/(2R). This 
flux is related to T by radiative diffusion, equation (J2J). The column depth is set 
by continuity to be y = Mt^/ (AitOslR 2 ) , where isL = h 2 jv is the timescale for 
viscous dissipation in the SL. The viscosity between the fresh material and the 
underlying WD is given by v « Th/(pv ( j > ) ot^v^h « osl^k^- Using an ideal 
gas equation of state, the above set of equations that can be solved for h. We 
set #sl ~ #disk using equation (j3J), since the disk may be setting the spreading 
angle as suggested by Figure 2. Equation (5) then gives a period 

P m =o - 60 s a disk 2 agL^Aj M 17 M 1 R 9 , (7J 

where Ai = A/10 -1 . Since this lowest order mode does not propagate in the 
azimuthal direction, we denote it with azimuthal wavenumber m = 0. Equation 
(7) provides both a scaling with M and a period suggestive of DNOs. 

The SL can contain additional modes, most notably those that propagate 
in the azimuthal direction. These modes have an observed frequency of w bs = 
\oj — muJsh\, where lo^l is the spin of the SL and m is the azimuthal wavenumber. 
Even though the layer is spinning quickly, Coriolis effects will not alter k, since 
6» SL £ w/(2u S l) (Bildsten, Ushomirsky & Cutler 1996). Setting P S L = 2vr/o;sL, 
where Psl ^, Pk, the next lowest order modes have 



Pm=±l 

for their periods. 



1 1 

P m =o T Ps[ 



(8) 



4.2. DNOs in the CV Population 

We compare the most recent compilations of DNO periods (Table 1 of Warner 
2004) and WD masses (Ritter & Kolb 2003) of CVs in Figure 3. We plot 
systems that have shown both high and low DNO periods as two separate points 
(these CVs are SS Cyg, CN Ori, and VW Hyi). The heavy, dashed line denotes 
the surface Keplerian period, Pk, for a given WD mass (using the mass-radius 
relation of Truran & Livio 1986). This demonstrates that DNOs exist both 
above and below so we consider SL nonradial oscillations with m = and 
m = —1. For each mode we plot periods for M = 10 16 — 10 18 g s _1 , the range 
expected during a DN outburst. To calculate the m = — 1 mode we must assume 
a period for the SL's spin, Psl ; so we plot mode periods for both Psl = Pk and 
Psl = 2Pk- This shows that our model is insensitive to the SL spin rate. 
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Figure 3. The DNO periods versus measured WD masses. The vertical 
bars show the range of periods observed from each CV, while the horizontal 
bars are mass measurement errors. CVs that have shows two distinct DNO 
periods are plotted as two separate points. The thick dashed line is the 
Keplerian period, The to — mode (wide, heavy shaded region) is 

plotted for M = 10 16 — 10 18 g s _1 (from top to bottom), using equation (7). 
The to = — 1 mode is also plotted for M = 10 16 — 10 18 g s -1 (from top to 
bottom), for both Psl = Pk (light shaded) and P$l = 2Pk (medium shaded), 
using equation (8). 

Besides predicting multiple classes of DNOs with different period ranges, 
we also predict that each of these groups should have a different dependence on 
M. This can be seen from the size of each shaded band, which is much wider 
for the m = mode than the m = — 1 mode. The DNOs with P ^ Pk show less 

variation with M, qualitatively consistent with the m = — 1 mode of our model. 
To further investigate the M dependence we plot our predicted DNO periods 
versus M for a M = I.OMq WD in Figure 4. This illustrates the shallower 
dependence for the m = — 1 mode. The separation and relative slopes of the 
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m = — 1 and m = modes are suggestive of the frequency doubling during the 
October 1996 outburst of SS Cyg (Mauche & Robinson 2001). 
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Figure 4. The SL nonradial oscillation period compared to M for a M = 
l.OM WD. The solid, heavy line is the m = mode, while the shaded region 
shows the m = — 1 mode for a range of Psl = (1-0 — 1.5) x Pk (bottom to 
top) . The dashed lines show the m = 1 mode for the same range of Psl • 



In Figure 4 we also plot the m = 1 retrograde mode for Psl = (1 — 1.5) x 
Pk- Due to the difference taken in equation (8), this mode has a complicated 
dependence on M, no longer being a power law nor monotonic. It has a higher 
period than typical DNOs, and may be relevant for the long period DNOS 
(lpDNOs; Warner, Woudt & Pretorius 2003). To identify these as modes requires 
checking for this predicted scalings between P and M. Such a test may negate 
a mode explanation for lpDNOs since Warner et al. (2003) claim to find no 
such correlation and thus favor a spin related mechanism. On the other hand, 
on separate occasions there have been 32 — 36 s (Robinson & Nather 1979) 
and 83 — 110 s (Mauche 2002b) oscillations seen from SS Cyg, in the domain 
expected for these modes. This wide spread of periods may be explained by the 
m = 1 mode's steeper dependence on M. The 32 — 36 s oscillations showed less 
coherence than typical DNOs, which could be due to these retrograde oscillations 
beating against the accretion disk. Oscillations from VW Hyi at ~ 90 s (Haefner, 
Schoembs & Vogt 1977) may be of similar origin. 
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An important property of WDs that would affect the modes is a magnetic 
field. A strong field inhibits shearing between the SL and WD and modifies the 
frequency of shallow surface waves. It is therefore interesting that no interme- 
diate polars (IPs) have shown DNOs or lpDNOs. Even LS Peg and V795 Her, 
neither of which are IPs, but both show polarization modulations (Rodriquez-Gil 
et al. 2001, 2002) indicative of a reasonably strong magnetic field, are without 
DNOs or lpDNOs. 

5. Discussion and Conclusions 

The SL model provides a new area of investigation in the study of accreting 
objects. The initial calculations by IS99 and in this paper describe the main 
features of spreading that will lead to further studies of how accreted material 
settles onto stars. This may include studying how differential rotation affects 
the underlying stellar surface. Spectral modeling of the SL, along comparisons 
with observations, will help in identifying if and when spreading is present. In 
this review we focus on comparisons with dwarf novae, but symbiotic binaries 
(Sokoloski 2004) and supersoft sources are also promising for studying the SL. 

We propose that DNOs in outbursting CVs are nonradial oscillations in the 
SL. A number of DNO properties are then simply understood: (1) the highly 
sinusoidal nature of the oscillations is consistent with nonradial oscillations, (2) 
the periods can change on the timescale of accretion because there is little mass 
in the layer (<; 10 21 g; PB04), (3) the periods vary inversely with M because they 
have the temperature scaling of shallow surface waves, (4) the covering fraction 
is naturally small for the SL, (5) the DNOs are only seen during DN outbursts 
when an optically thick layer of material builds up at the equator, and (6) the 
largest pulsed amplitude is in the EUV, consistent with the SL temperature. 

There are difficulties that must be answered about this idea for explaining 
DNOs. From the SL model we borrowed the concept of hot material in hy- 
drostatic balance covering a small fractional area of the WD, but we ignored 
important details of these calculations, and instead only presented a very simple 
model for the modes. The data do not require such additional complications, so 
we refrain from including them for now. In a more sophisticated model it would 
probably still be true that k oc 1/(9slR), but an eigenvalue calculation would 
determine the constant of proportionality. 

Further studies should work toward understanding the excitation mecha- 
nism for the modes, which would explain the high coherence typical of DNOs. 
Material deposited at the WD equator spreads quickly, ~ 10 — 100 s (PB04a), so 
that following an outburst, when M has fallen, the accreted material will spread 
over the star and not be seen as a separate hot component. This also means 
that the material is moving through the oscillating region on timescales of order 
the mode period. This short timescale is problematic if the modes are to remain 
coherent, which puts limits any proposed excitation mechanism. 

Our explanation of DNOs raises interesting questions about the relation- 
ship between oscillations originating from accreting compact objects. Mauche 
(2002b) showed that there is a correlation in the high to low oscillation frequency 
ratio of WDs, NSs, and black holes (BH). Using this picture, DNOs are asso- 
ciated with the kilohertz QPOs of low mass X-ray binaries. Suggestively, the 
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Fourier frequency resolved spectroscopy of NSs (Gilfanov, Revnivtsev &: Molkov 
2003) imply that both the normal branch oscillations and the kilohertz QPOs 
originate in the NS BL. Using a NS mass and radius our model provides fre- 
quencies in the range expected for kHz QPOs, but it does not explain their P-L 
relation (the "parallel tracks"; van der Klis 2000). Interestingly, DNOs may 
also show the parallel tracks phenomenon, as seen in three observations of SS 
Cyg (Mauche 2002a), once again supporting the correlation. On the other hand, 
continuing the analogy to WD and BH oscillations is problematic because in the 
case of BHs there is no surface for nonradial oscillations. 
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